X 

5-! 



Stellar Clusters & Associations: A RIA Workshop on Gaia 

Granada, 23 rd - 2f h May 2011 



Feedback Regulated Star Formation: From Star 
Clusters to Galaxies 



Sami Dib 1 

1 Astrophysics Group, Blackett Laboratory, Imperial College London, London, SW7 2AZ, 
United Kingdom; s.dib@imperial.ac.uk 



Abstract 



O 
(N 

3 

in 

< 
O 

This paper summarises results from semi-analytical modelling of star formation in proto- 
cluster clumps of different mctallicitics. In this model, gravitationally bound cores form 

+3 uniformly in the clump following a prescribed core formation efficiency per unit time. After 

a contraction timescale which is equal to a few times their free-fall times, the cores collapse 

i__i into stars and populate the IMF. Feedback from the newly formed OB stars is taken into 

account in the form of stellar winds. When the ratio of the effective wind energy of the 
winds to the gravitational energy of the system reaches unity, gas is removed from the clump 

vq and core and star formation are quenched. The power of the radiation driven winds has 

a strong dependence on metallicity and increases with increasing metallicity. Thus, winds 
from stars in the high metallicity models lead to a rapid evacuation of the gas from the 
protocluster clump and to a reduced star formation efficiency, SFE exp , as compared to 

£■ — their low metallicity counterparts. By combining SFE exp with the timescales on which gas 

expulsion occurs, we derive the metallicity dependent star formation rate per unit time in 
this model as a function of the gas surface density S ff . This is combined with the molecular 
gas fraction in order to derive the dependence of the surface density of star formation T,sfr 
on X g . This feedback regulated model of star formation reproduces very well the observed 
star formation laws extending from low gas surface densities up to the starburst regime. 
Furthermore, the results show a dependence of "Ssfr on metallicity over the entire range 



of gas surface densities, and can also explain part of the scatter in the observations. 

1 THE SFE IN STELLAR CLUSTERS AND IN GALAXIES 

The star formation efficiency (SFE) is one of the essential quantities that regulates the dy- 
namical evolution and chemical enrichment of stellar clusters, the interstellar medium, and 
galaxies (e.g., Boissier et al. 2001; Krumholz & McKee 2005; Dib et al. 2006; Dib et al. 
2009; Dib et al. 2011). The SFE is commonly defined as being the fraction of gas which is 
converted into stars in a system of a given mass, be it a protocluster molecular clump, an 
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Figure 1: Global galactic star formation efficiency SFE ga i as a function of the galactic global 
metallicity. The data for M33 and NGC 6822 are from Gratier et al. (2010a,b), for the SMC from 
Leroy et al. (2006), and the average spiral value from the sample of Murgia et al. (2002). Adapted 
from Dib ct al. (2011) 

entire giant molecular cloud (GMC), or a galaxy. In an isolated clump, the time evolving 
SFE is usually defined as being: 



SFE{t) 



M duster (t) 



M. 



clump 



(1) 



where M c i U ster(t) is the cluster mass at time t and M c i U mp is the initial clump mass. In 
observed protocluster clumps with embedded or semi-embedded clusters, the SFE is usually 
approximated by: 



SFE, 



obti 



^cluster 
Mduster + M gas 



(2) 



where M gas is the mass of the star forming gas. Measured SFEs of nearby embedded clusters, 
using Eq.[2]yield values that fall in the range 0.1-0.5 (Wilking & Lada 1983; Rengarajan 1984; 
Wolf et al. 1990; Pandey et al. 1990; Lada et al. 1991a,1991b; Warin et al. 1996; Olmi & 
Testi 2002). These SFE values are larger than those obtained for entire GMCs which are 
observed to fall in the range of a few percent (e.g., Duerr et al. 1982, Fukui & Mizuno 1991; 
Evans et al. 2009). On galactic scales, the star formation efficiency is usually defined as 
being the inverse of the molecular gas consumption time and is given by: 
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Figure 2: Time evolution of the fraction of cloud mass of the cloud that is in dense cores in two 3D 
numerical simulations of magnetized, turbulent, and self gravitating isothermal molecular cloud (Dib 
et al. 2010a). One simulation has a mass-to magnetic flux ratio of \ib — 2.8 (mildly supercritical) 
and the other has hb = 8.8 (strongly supercritical). The upper corresponds to cores detected using a 
density threshold of n t f lr — 20 n aver and the lower panel corresponds to cores detected at the density 
threshold of 160 n aver , where n aver — 500 cm~ 3 is the average number density in the cloud. The cores 
identified at the threshold density of 160 at , er corresponds to a population of gravitationally bound 
cores. Adapted from Dib et al. (2010a) 
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where SFR ga i is the galactic star formation rate (in M yr 1 ), and M# 2 the mass of the 
molecular hydrogen gas. In large spiral galaxies SFE ga i ~ 0.5 x 10~ 9 yr -1 (e.g., Kenni- 
cutt 1998, Murgia et al. 2002; Leroy et al. 2008). Fig. [I] displays the global SFE ga i for a 
few selected nearby galaxies as a function of their global metallicity. The trend in Fig. [T] is 
strongly suggestive of an SFE ga i which increases with decreasing metallicity. Several phys- 
ical processes regulate the core formation efficiency (CFE) within a star forming molecular 
clump/cloud, and as a consequence, the SFE. Supersonic turbulence can provide support 
against the global collapse of the clump/cloud, or at least delay it (e.g., Vazquez-Semadeni 
& Passot 1999). However, on scales smaller than the energy injection scale, but larger than 
the sonic scales in the cloud, supersonic turbulence produces local compressions (i.e., cores) 



Feedback Regulated Star Formation 



o 

GO 



"O 



~TD 



o 








40 



60 



80 



ivi a s s 

Figure 3: Stellar mass loss rates for stars in the mass range 5-80 M on the main sequence, and 
for various metallicities. The stellar mass loss rates have been calculated using the stellar character- 
istics (effective temperature, stellar luminosity and radius) computed using the stellar evolution code 
CESAM coupled to the stellar atmosphere model of Vink et al. (2001). Over-plotted to the data are 
fourth order polynomial fits. The parameters of the fit functions can be found in Dib et al. (2011). 
Adapted from Dib et al. (2011) 



(e.g., Padoan 1995; Vazquez-Semadeni et al. 2003; Dib et al. 2010a) of which a fraction can 
be 'captured' by gravity and proceed to collapse into stars (e.g., Vazquez-Semadeni et al. 
2005a; Dib et al. 2007a; Dib & Kim 2007; Dib et al. 2008a; Dib et al. 2010a). Krumholz 
&: McKee (2005) formulated an analytical theory in which the core formation efficiency per 
free-fall time, CFEff, is shown to decrease with an increasing sonic Mach number and an 
increasing virial parameter and is given by CFEjf « 0.15a~ r ' .M -0 ' 32 . 

Magnetic fields play an important role in determining the fraction of gravitationally 
bound gas in star forming clouds/clumps. Results from numerical simulations show that 
stronger magnetic fields (in terms of magnetic criticality) lower the rate of dense core forma- 
tion in a star forming molecular clump/cloud (e.g., Vazquez-Semadeni et al. 2005b; Price & 
Bate 2008; Dib et al. 2008a; Li et al. 2010; Dib et al. 2010a). Dib et al. (2010a) showed that 
the CFE per unit of the free-fall time of the cloud, CFEff, are of the order of ~ 6 % and 
~ 33 % for clouds with mass to-magnetic flux ratios of \x = 2.2 and 8.8, respectively (with 
[i being normalised by the critical mass-to-flux ratio for collapse, see Fig. ^1). The role of 
stellar feedback in setting the final value of the SFE has been investigated by several authors. 
The role of protostellar outflows has been studied theoretically by Adams & Fatuzzo (1996), 
Matzner & McKee (2000) and numerically by Nakamura & Li (2007) and Li et al. (2010). 
Supernova explosions are an efficient way of removing gas from the protocluster region (e.g., 
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Figure 4: The power of the stellar winds, or wind luminosities, for stars in the mass range 5-80 M Q 
on the main sequence, and for various metallicities. The stellar mass loss rates have been calculated 
using the stellar characteristics (effective temperature, stellar luminosity and radius) computed using 
the stellar evolution code CESAM coupled to the stellar atmosphere model of Vink et al. (2001). The 
values of Voo have been calculated using the derivation by Leitherer et al. (1992). Over-plotted to 
the data are fourth order polynomials. The parameters of the fit functions can be found in Dib et al. 
(2011). Adapted from Dib et al. (2011) 



Parmentier et al. 2008; Baumgardt et al. 2008). However, they occur after a few million 
years from the time massive stars have formed. Another form of stellar feedback is associated 
with O and B stars, in their main sequence phase, and eventually beyond. OB stars emit UV 
radiation which ionises the surrounding gas and heats it to temperatures of ~ 7000 — 10 4 K. 
This warm and ionised bubble provides the environment in which particles accelerated from 
the stellar surface by interaction with some of the stellar radiation propagate outwards. In 
the following sections, we present a model which describes the co-evolution of the mass func- 
tion of gravitationally bound cores and the IMF in a protocluster region. In the model, dense 
cores form in the protocluster clump uniformly in time following a specified core formation 
efficiency per unit free-fall time of the clump, CFEff. The dense cores have lifetimes of a 
few times their free-fall times, after which they collapse to form stars and populate the IMF. 
Stellar winds from the newly formed massive stars (M* > 5 M ) inject energy into the clump 
and when the ratio of the effective wind energy to the gravitational energy of the clump 
reaches unity, gas is removed from the clump and core and star formation are quenched. We 
discuss the dependence of the final star formation efficiency on the metallicity. Finally, we 
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present a derivation of the star formation laws in galaxies in the context of this metallicity 
dependent, feedback regulated model of star formation. 

2 The Model 

2.1 Protocluster Clumps 

Several studies have established that star clusters form in dense (> 10 3 cm~ 3 ) clumps em- 
bedded in a lower density parental molecular cloud (e.g., Lada & Lada 2003; Klein et al. 
2005; Rathborne et al. 2006; Kauffmann & Pillai 2010; Csengeri et al. 2011; Hernandez et 
al. 2011). Saito et al. (2007) studied, using the C 18 molecular emission line, a large sample 
of cluster forming clumps whose masses and radii vary between [15-1500] M and [0.14-0.61] 
pc, respectively. The mass-size and velocity dispersion-size relations of clumps in the sample 
of Saito et al. (2007) were fitted by Dib et al. (2010b), obtaining: 

M dump (M Q ) = 10 3 - 62±014 J R c 2 - 54±0 - 25 (pc), (4) 

and 

t> c (km s- 1 ) = 10 - 45±0 - 08 i?°- 44±0 - 14 (pc), (5) 

where M c i ump is the mass of the clump, R c its radius, and v c the scale dependent gas velocity 
dispersion. In this work, we adopt a protocluster clump model that follows an r~ 2 density 
profile: 

Pc{r) = rrwkop (6) 

where R c q is the clump's core radius, p c o is the density at the centre. For a given value of 
Mciump, the central density p C Q is given by: 

Pc ° ~ AnRl [{R c /R c0 ) - arctan(i? c / J R c0 )] ' [ ' 

The temperatures of the cluster forming clumps are observed to vary between 15 and 
70 K (e.g., Saito et al. 2007; Rathborne et al. 2010). In order to further constrain the models 
and minimise the number of parameters, we relate the sizes of the protocluster clumps to 
their masses using the mass-size relation of Saito et al. In the absence of detailed information 
about the velocity dispersion inside the cores in the Saito et al. (2007) study, we assume that 
the clump-clump velocity dispersion they derived (i.e., Eqrc3j) is also valid on the scale of the 
clumps themselves and of their substructure. 

2.2 Prestellar Cores Model 

Whitworth &; Ward-Thompson (2001) applied a family of Plummer sphere-like models to 
the contracting prestellar dense core L1554, which is representative of the population of 
gravitationally bound cores in clumps that are considered in this work. They found a good 
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agreement with the observations of L1554 if the density profile of the core has the following 
form: 

Pp{rp) = [l + (r p /i? p o) 2 ] 2 ' (8) 

where p p o and R p q are the central density and core radius of the core, respectively. Note that 
the radius of the core, R p , depends both on its mass and on its position within the clump. 
The dependence of R p on r requires that the density at the edges of the core equals the 
ambient clump density, i.e., p p {R p ) = p c ( r )- This would result in smaller radii for cores of 
a given mass when they are located in their inner parts of the clump. The density contrast 
between the centre of the core and its edge is given by: 



C(r) = -^ = M 

Pc{r) Pco 



r 
-KcO 



(9) 



Depending on its position r in the clump, the radius of the core of mass M, R p , can 
be calculated as being R p (r,M) = a(r) R p o(r,M), where: 

( M \ 1/3 / air) \~ 1/3 

^(r,M)=^— j (arctan[a(r)] - 1 ^ F ) , (10) 

and with a(r) = (C^r) 1 ' 2 — l) 1 ' 2 . With our set of parameters, the quantity C 1 ' 2 — 1 is 
always guaranteed to be positive. The value R p (r, M) can be considered as being the ra- 
dius of the core at the moment of its formation. The radius of the core will decrease as 
time advances due to gravitational contraction. We assume that the cores contract on a 
timescale, t c <mt,p which we take to be a few times their free fall timescale tff iCore , and which 

1 /9 

is parametrized by t con t, p (r, M) = v tff tCore (r,M) = iy(3n/32 Gp~ p (r, M)) ' , where G is the 
gravitational constant, v is a constant > 1 and p~ p is the radially averaged density of the 
core of mass M, located at position r in the clump. The time evolution of the radius of 
a core of mass M, located at position r in the cloud is given by a simple contraction law 
Rp(r,M,t) = Rp(r,M,0) e ~('/ tco ™ t -p). Both observational (Jessop & Ward-Thompson 2000; 
Ward-Thompson et al. 2007) and numerical (Vazquez-Semadeni et al. 2005a; Galvan-Madrid 
et al. 2007; Dib et al. 2008c; Gong & Ostriker 2011) estimates of gravitationally bound cores 
lifetimes tend to show that they are of the order of a few times their free-fall time. One 
important issue is the choice of the cores central density, p P Q. In this work, we first assume 
that the minimum density contrast that exists between the centre of the core and its edge is 
of the order of the critical Bonnor-Ebert value and that is > 15. Secondly, we assume that 
the density contrast between the centre and the edge of the cores depends on their masses 
following a relation of the type p p o oc M M . Thus, the density contrast between the centre and 
the edge for a core with the minimum mass we are considering, M m i n (typically M m i n = 0.1 
M ) is 15, whereas for a more massive core of mass M, the density contrast will be equal to 
15 x (M/M m i n )^. Observations show that p varies in the range [0 — 0.6] (e.g., Johnstone & 
Bally 2006). 
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Figure 5: Time evolution of the pre-stellar core mass function (left), and stellar mass function (right) 
in the protocluster clump with the fiducial model parameters. The last time-step shown is 1.79 £//, c 
which corresponds to the epoch at which gas is expelled from the protocluster clump. Adapted from 
Dib ct al. (2011) 

2.3 The Initial Prestellar Cores Distributions 

We assume that dense cores form in the clump as a result of its gravo-turbulent fragmenta- 
tion. The mass distributions of the core formed at every epoch is described using the formu- 
lation of Padoan & Nordlund (2002). Thus, the local mass distributions of cores, N(r,M), 
are given by: 



N(r,M) dlog M = f (r) M -3 ^ 4- ^ 

-M 

P(Mj)dMj dlog M, 



("J 



where /3 is the exponent of the kinetic energy power spectrum, E^ oc k &, and is related to 
the exponent a of the size- velocity dispersion relation in the clump with f3 = 2a + 1. The 
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Figure 6: Time evolution of the ratio of the effective kinetic energy generated by stellar winds and 
the gravitational energy of the protocluster clump. Time is shown in units of the protocluster clump 
free fall timescale £//. c . The horizontal dashed corresponds to Ek^ W ind/E grav = 1 with k = 0.1. The 
full line correspond to the fiducial model with the CFEff = 0.2 and the dashed and triple dot-dashed 
to cases with CFEff =0.1 and CFEff — 0.3, respectively. Diamonds correspond to the epochs at 
which the gas is evacuated from the cluster in the three models and the processes of core and star 
formation arc terminated. Adapted from Dib et al. (2011). 

local normalisation coefficient fo(r) is obtained by requiring that J M max N(r,M) dM = 1 
in a shell of width dr, located at distance r from the clump's centre. P{Mj) is the local 
distribution of Jeans masses given by: 



P{Mj) dMj = 2 -^hMf 



2naj 



exp 



1 /In Mj- 

n o7 



A 



dM, 



(12) 



where Mjq is the Jeans mass at the mean local density, and a^ is the standard deviation of 
the density distribution which is a function of the local thermal rms Mach number. Therefore, 
the local distribution of cores generated in the clump, at an epoch r, N(r, M, r), is obtained 
by multiplying the local normalised function N(r, M) by the local rate of fragmentation such 
that: 



N(r, M, T )dt 



CFE ff {r)p c {r) 



dt 



<M>{r) t cont)P (r,M)t ffy 



N(r,M), 



(13) 



where dt is the time interval between two consecutive epochs, < M > is the average core 
mass in the local distribution and is calculated by < M >= J^ max M N(r, M, 0) dM, and 
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CFEff is a parameter smaller than unity which describes the local mass fraction of gas that 
is transformed into cores per free fall time of the protocluster clump, t//, c - In the present 
study, we assume that CFEff is independent of r. 
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Figure 7: Time evolution of the SFE in the protocluster clump. The full line corresponds to the 
fiducial model with the CFEff = 0.2 and the dashed and triple dot-dashed to cases with CFEff = 0.1 
and CFEff = 0.3, respectively. Time is shown in units of the protocluster clump free fall timescale 
i//, c - Diamonds correspond to the epochs at which the gas is evacuated from the cluster and indicates 
the final value of the SFE, SFEf, in the three models. Adapted from Dib et al. (2011). 

In this model, the formation of cores in the protocluster clump, and consequently star 
formation, are terminated whenever the fraction of the wind energy stored into motions that 
oppose gravity exceeds the gravitational energy of the clump. Thus, at any epoch t < t exp 
(t exp is the epoch at which gas is expelled from the clump), gas is removed from the clump 
only to be turned into stars. We take into account the feedback generated by the stellar 
winds of massive stars (M* > 5 Mq). In order to calculate reliable estimates of the feedback 
generated by metallicity dependent stellar winds, we proceed in two steps. In the first step, 
we use a modified version of the stellar evolution code CESAM (see appendix 1 in Piau et al. 
2011) to calculate a grid of main sequence stellar models for stars in the mass range [5-80] M 
(with steps of 5 M ) at various metallicities Z/Z & = [1/10, 1/6, 1/3, 1/2, 1, 2] (Z & = 0.0138). 
The evolution of massive stars is followed using the CESAM code for ~ 1 Myr , on the main 
sequence. The characteristic stellar properties, which are the effective temperature T e ff, the 
luminosity L*, and the stellar radius R+ are then used in the stellar atmosphere model of 
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Figure 8: Time evolution of the ratio of the kinetic energy generated by stellar winds and the 
gravitational energy of the protocluster clump. Time is shown in units of the protocluster clump free 
fall timescale i//, c - The horizontal dashed corresponds to Ek,wind/E grav = 1 with k = 0.1. Diamonds 
correspond to the epochs at which the gas is evacuated from the cluster in the three models. The 
mass of the clump in these models is 10 5 M Q . Adapted from Dib et al. (2011). 

Vink et al. (2001) in order to calculate the stellar mass loss rate M*. Vink et al. (2001) did 
not derive the values of the terminal velocities of the winds (i><x>), therefore, we use instead 
the derivations of Voo obtained by Leitherer et al. (1992). Fig. [3] displays the mass loss rates 
calculated for OB stars at the various metallicities. The power of the stellar winds is given 
by M+v^. This quantity is displayed in Fig. UJ for the models with different metallicities. 
Both log(M*) and M*t^ are fitted with fourth order polynomials (overplotted to the data) 
and whose coefficients are provided in Dib et al. (2011). The M+v^ — M* relations displayed 
in Fig. [4] allow for the calculation of the total wind energy deposited by stellar winds. The 
total energy from the winds is given by: 



Jt'=0 JM t =5 M n \ 2 / 



(14) 



We assume that only a fraction of E W i n d will be transformed into systemic motions that 
will oppose gravity and participate in the evacuation of the bulk of the gas from the proto- 
cluster clump. The rest of the energy is assumed to be dissipated in wind-wind collisions or 
escape the wind bubble. The effective kinetic wind energy is thus given by: 



^k.wind — ^ J^wind: 



(15) 
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where k is a quantity < 1. It is currently difficult to estimate k as its exact value will vary 
from system to system depending on the number of massive stars, their locations, and their 
wind interactions. As a conservative guess for the fiducial model, we take k = 0.1. E^^wind is 
compared at every timestep to the absolute value of the gravitational energy, E grav , which is 
calculated as being: 

16 r Rc 

Egrav = -^tt 2 G \ p c (r) 2 r A dr, (16) 



where p c is given by Eq. |6j 
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3.1 The co-evolution of the core mass function and of the IMF 

Whenever a population of cores of mass M, located at a distance r from the centre of the 
clump has evolved for a time that is equal to its contraction timescale, it is collapsed into 
stars. Thus, the local number of cores of a given mass, at a given epoch r, is the sum of all 
the local populations of cores of the same mass that have formed at all epochs anterior or 
equal to the considered epoch with the additional step of subtracting from that sum the cores 
of the same mass that have readily collapsed into stars. The local populations of cores of 
various ages are evolved separately as they are each in a different phase of their contraction, 
and will collapse and form stars at various epochs. Thus, the total local number of cores of 
a given mass M, at a time t, is given by: 

N(r,M,t) = J2 N (r,M, n ,t). (17) 

n<t 

We assume that only a fraction of the mass of a core ends up locked in the star. We 
account for this mass loss in a purely phenomenological way by assuming that the mass of a 
star which is formed out of a core of mass M is given by M* = (,M, where £ < 1. Matzner & 
McKee (2000) showed that £ can vary in the range 0.25 — 0.75. In this section, we describe the 
results for a fiducial model. In this model, the mass of the clump is taken to be M c ; ump = 10 5 
M , the metallicity is Z = Z@, and The CFE is CFEff = 0.2, where tffd is the free-fall time 
of the clump and is given by tff yC = (37r/32 Gpc) 1 ' 2 , and p c the average density of the clump. 
The other quantities have been taken to be equal to the most commonly cited observational 
determinations and have been set, in the fiducial model as well as in all other models, to the 
following values: R c o = 0.2 pc, v = 3, p = 0.2, k = 0.1. We also assume that £ = 1/3. Fig. [5] 
displays the time evolution of the CMF (left column) and of the IMF (right column) in the 
fiducial model. By t ~ 0.3 tft c , the first stars form. In this model, since p = 0.2 > 0, the 
most massive stars form first, as the most massive cores tend to be, at any given position 
in the clump, more centrally peaked and thus have shorter lifetimes. As time advances, the 
IMF becomes fully populated. On the other hand, the CMF ceases to evolve (i.e., there is 
no accretion or coalescence in this model as in Dib 2007, Dib et al. 2007b; Dib et al. 2008b; 
and Dib et al. 2010a) as the numbers of cores that are newly formed at each position in 
the protocluster clump is balanced by an equal number of cores which collapses and forms 
stars. The final IMF of the cluster is established at ~ 1.79 tff c , which corresponds to the 
last epoch shown in Fig. [5] At this epoch, Fig. M shows that the ratio Ek^ wina \/ E grav reaches 
unity (full line), and as a consequence, gas is expelled from the protocluster region. Fig. u\ 
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displays the time evolution of the SFE in this model (full line). At t = t exp = 1.79 tff )C , the 
final value of the SFE is SFE exp ~ 9.05 x 1(T 2 . 

In this fiducial model, we have adopted a CFE value of CFEff = 0.2. This is an 
intermediate value between the values of ~ 0.06 and ~ 0.33 measured by Dib et al. (2010a) 
in numerical simulations of molecular clouds with two different degree of magnetisation. At 
first glance, it may appear that increasing the CFEff by a given factor will lead to an 
increase in SFE exp by approximately the same factor. However, for a fully sampled IMF, 
a larger CFEff value implies that a larger number of OB stars will be formed and deposit 
larger amounts of feedback by stellar winds in the protocluster region. This in turn leads to 
a faster evacuation of the gas and to a limitation of the SFE exp . Fig. [6] shows that in two 
other models similar to the fiducial case but with a different CFEff, the evacuation of the 
gas occurs faster for the higher CFE case (i.e., case with CFEff = 0.3) and slower for the 
low CFE case (i.e., CFEff = 0.1). The time evolution of the SFE in these two additional 
models is compared to the fiducial case in Fig. [7| The final values of the SFE at t = t exp 
are SFE exp ~ 6.59 x 10~ 2 , ~ 9.05 x 10~ 2 , and ~ 0.11 for the cases with CFEff = 0.1,0.2, 
and 0.3 respectively. A factor of 2 and 3 increase in the CFEff leads to an increase of the 
SFE exp only by factors of ~ 1.37 and ~ 1.64, respectively. This implies a strong regulation 
of the effect of a varying CFE by stellar feedback on the resulting SFE exp in protocluster 
clumps. 

3.2 The dependence of the star formation efficiency on metallicity 
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Figure 9: S'i 7 '-E ea ;p-Metallicity relations for clumps with the fiducial mass of 10 5 M and for different 
core formation efficiencies of CFEff = 0.1, 0.2, and 0.3. Over-plotted to the data are fit functions 



(Eq. 25) whose parameters are given in Dib et al. (2011). Adapted from Dib et al. (2011) 
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Figure 10: Dependence of the quantities SFE exp (final star formation efficiency) and n, 

(ratio of the expulsion time to the free-fall time) for selected values of the protocluster clump masses 

and metallicies. Adapted from Dib (2011). 

Fig. [8| displays the time evolution of the ratio Ef.,wind/ ^grav in models similar to the 
fiducial case but with metallicities varying between Z = 0.1 Zq and Z = 2 Zq. All six models 
have the fiducial value of the CFE, CFEft = 0.2. In models with lower metallicities, the 
power of the stellar winds in weaker (i.e., Fig. El) and the evacuation of the gas occurs at later 
epochs as compared to the higher metallicity cases. In the model with Z = 2 Zq, the gas is 
expelled from the protocluster region at t ~ 1.4 tft c while in the model with Z = 0.1 Zq, 
the gas expulsion is delayed until t ~ 4.5 i//. c - For a given CFEff, longer timescales imply 
a larger SFE exp . These values are plotted versus metallicity in Fig. [9| A clear trend is 
observed in which the SFE exp increases with decreasing metallicity (diamonds). We fit the 
SFE exp — Z points in Fig. [9] with the following functional form: 



SFE f = C z e 



^ lo <4) 



(18) 



Fo the case case with CFEfj=0.2, the fit parameters are Cz = 0.091 ± 6.8 x 10 -4 and 
tz = 0.91 ± 7.9 x 10~ 3 . We have also repeated the calculations at the various metallicities 
using the additional values of CFEff = 0.1 and CFEff = 0.3. The dependence of the 
SFE exp on metallicity in these cases are also displayed in Fig. [9| Fig. M shows that the 
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SFE exp varies typically by a factor of ~ 1.3 — 1.4 and ~ 1.6 — 1.7 for variations of the CFE 
by factors of 2 and 3, respectively. This argues for a strong regulation of the SFE by feedback. 

4 The Kennicutt- Schmidt laws in the feedback regulated model 
of star formation 

Using the above described model, it is possible to derive a metallicity dependent star forma- 
tion rate (SFR) at a given surface density T, g . For a summary on the recent observational 
determinations of the star formation laws see Kennicutt (2008), Elmegreen (2011), and Dib 
(2011). The star formation rate surface density in the feedback regulated mode of star for- 
mation is given by: 



(SFE exp 
\texp 



^SFR = S 3 /ff2 77 v j (19) 



where (SFE exp ) and (t exp ) are, respectively, the characteristic SFE exp and the epoch at 
which gas is expelled from the protocluster region for the clump mass distribution associated 
with a given S ff . Writing (t exp ) in terms of the clumps free-fall time (£//), Eq. M becomes: 



^sfr = S 9 fn 2 



(SFE exp ) 



(nexp) (iff) 



Z g fH 2 { -^- (20) 



where /*,// is the dimensionless star formation efficiency and which corresponds to the mass 
fraction of the molecular gas that is converted into stars per free-fall time i// jC of the clumps 
and. (/*,//) and (£//) represent characteristic values of /*,// and tff )C for the spectrum of 
clump masses found in the giant molecular for a given value of E g . The quantity fn 2 is the 
mass fraction of the total gas that is in molecular form. In this work, we use the functional 
form of fjj 2 obtained by Krumholz et al. (2009a) who derived fn 2 as a function of the 
gas surface density and metallicity (see their paper or Dib 2011 for the detailed formula). 
{tffc} can be approximated by the free-fall time of the clump with the characteristic mass 

tff, c (M char ) = 8S^ 3/ M c [ arfi Myr where M charfi = M char /10 6 M . The characteristic mass 
M char is given by : 

M char = / M d N(M cl )dM ch (21) 

" ^" cl ,min 

where N(M c i) is the mass function of protocluster forming clumps which we take to be 
N(M c i) = A c iM^ , and Ad is a normalisation coefficient given by A c i J M . cl ' max ' GMC N(M c {)dM c i 
e, where < e < 1 is the mass fraction of the GMCs that is in protocluster clumps at any 
given time. In this work we use e = 0.5. The minimum clump mass M^ TO j„ is taken to be 
2.5 x 10 3 M (this guarantees, for final SFEs in the range of 0.05-0.3 a minimum mass for 
the stellar cluster of ~ 50 M ) and the maximum clump mass is 10 8 M Q . The characteristic 
GMC mass is determined by the local Jeans mass and is given by: 
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Figure 11: Star formation efficiency per unit free-fall time in the protocluster clump in the metallicity- 
dependent feedback model. The star formation efficiencies per free-fell use a core-to-star efficiency 
conversion factor of 1/3. The left panel displays /*,// as a function of both M c i and Z = Z/Zq in 
the original data, while the right panel displays the analytical fit function to this data set given in 
Eq.|8] Adapted from Dib (2011). 



M GM c = 37 x 10 b 



85 M pc" 



M<7 



©• 



(22) 



Fig. 12 (top) displays M c h ar as a function of S g . Fig. 10 displays the dependence of 



SFE exp and t exp on clump mass and metallicity for a series of models in which the CFEff = 
0.2. The SFE exp depends strongly on metallicity but weakly on mass whereas t exp displays 
a clear dependence on both quantities. The quantity /*,// = SFE exp /n exp is displayed in 
Fig. [ill (left panel) as a function of mass and metallicity (Z = Z/Zq). A fit to the /*,// 
(Md, Z ) data points with a 2-variables second order polynomial yields the following relation 



(shown in Fig. 11, right panel): 



Ujf(M d , Z ) = 11.31 - 4.311og(M cl ) + 0.41[log(Mci)] 

-8.28Z + 3.20Z'log(M cl ) - 0.32Z'[log(M c i)] 

+2.30Z' 2 - 0.89Z' 2 log(M cl ) + 0.08Z' 2 [log(M cl )] :; 



Using Eq. 23, it is then possible to calculate (/*,//): 

rmax(M clmax ,M G Mc) 



(f*,ff)(Z,X 9 ) 



M r 



h Jf (M ch Z)N(M cl )dM d . 



(23) 



(24) 



Fig. 12 (bottom) displays (/*,//) (Z , S fl ) for values of Z in the range [0.1 — 2]. We 
assume that there is a critical value of E 9 = 85 Mq pc~ 2 below which clumps are pressurised 
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Figure 12: Characteristic clump mass as a function of the gas surface density (Eq. Ul top panel) and 
the star formation efficiency per unit free-fall time in this feedback regulated model of star formation. 
Adapted from Dib (2011). 

by their internal stellar feedback, such that T, c i = T, g ^ cr it where T, g < T, g>cr it and X c ; = 
^GMC = ^ 9 when E 9 > £ SjCr it- With the above elements, the star formation law can be 
re- written as: 



^SFR = Y^fH 2 {^g,c,Z )H g 



ML K 

{f*jf)(z') ( - 



• 9 <r 1 

' 85 M pc" 2 ^ ^ 



M 



1/4 
char, 6 



( ^ \ 

V85 M Q pc- 2 J 



3/4 



(25) 



' 85 M pc- 2 



> 1 



where ^sfr is in M Q yr x 
24 



kpc 2 , M char is given by Eq. 21, and (/*,//) by Eqs. 



23 



and 
Fig. 13 (top panel) displays the results obtained using Eq. 25 for S 9 values starting 
from low gas surface densities up to the starburst regime. The results are calculated for the 
metallicity values of Z = [0.1, 0.3, 0.5, 1, 2]. The results are compared to the sub-kpc data of 
Bigiel et al. (2008,2010) and to the normal and starburst galaxies results of Kennicutt (1998) 
(Fig. 13, bottom panel displays the results of Krumholz et al 2009 which is based on the 
regulation of star formation by turbulence in GMCs). Our models fits remarkably well the 
observational results over the entire range of surface densities. Furthermore, the segregation 
by metallicity extends beyond the low surface density regime up to the starburst regime 
where a segregation in metallicity of ~ 0.5 dex is observed, in contrast to the Krumholz et al. 
(2009b) models which do not contain a metallicity dependence in the intermediate to high 
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Figure 13: Star formation laws in the feedback-regulated star formation model (this work, top panel), 
and in the Krumholz et al. (2009a) model (bottom panel). Overplotted to the models are the normal 
and starburst galaxies data of Kennicutt (1998) and the combined sub-kpc data (4478 subregions in 
total) for 11 nearby galaxies from Bigiel et al. (2008,2010). The Bigiel et al. (2008,2010) data is shown 
in the form of a 2D histogram with the color coding corresponding, from the lighter to the darker 
colors to the 1,5,10,20, and 30 contour levels. The displayed theoretical models cover the metallicity 
range Z = Z/Z Q = [0.1,2]. Adapted from Dib (2011). 

surface density regimes. Furthermore the solar metallicity curve in our model overlaps with 
a significant fraction of the sub-regions in the data of Bigiel et al. (2008,2010) in contrast to 
the Krumholz et al. model. 



5 CONCLUSIONS 



We have described a model for star formation in protocluster clumps of difference metallicities. 
The model describes the co-evolution of the dense core mass function and of the IMF in the 
clumps. Cores form uniformly over time in the clumps following a prescribed core formation 
efficiency per unit time. Cores contract over timescales which are a few times their free fall 
time before they collapse to form stars. Feedback from the newly formed OB stars (> 5 M ) 
is taken into account and when the ratio of the cumulated effective kinetic energy of the winds 
to the gravitational energy of the system (left over gas+stars) reaches unity, gas is expelled 
from the clump and further core and star formation are quenched. The radiation driven winds 
of OB stars are metallicity dependent. Metal rich OB stars inject larger amount of energy 
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into the clump than their low metallictiy counterparts and thus help expel the gas on shorter 
timescales. This results in reduced final star formation efficiencies in metal rich clumps in 
comparison to their low metallicity counterparts. Both the final star formation efficiency and 
the gas expulsion timescales are combined for a grid of clump models with different masses 
and metallicities in order to calculate the star formation efficiency per unit time (/*,//) in this 
feedback regulated model of star formation. We calculate the characteristic value of /*,// for 
a clump mass distribution associated with a gas surface density, T, g . This is combined with a 
description of the molecular mass fraction as a function of E 9 and the assumption that there 
is a critical surface gas density (£ 9 = 85 M pc~ 2 ) above which the protocluster clumps 
and their parent giant molecular clouds switch from being pressurised from within by stellar 
feedback to being confined by the external interstellar medium pressure. The combination 
of these three elements allows us to construct the star formation laws in galaxies going from 
low gas surface densities up to the starburst regime. Our models exhibit a dependence on 
metallicity over the entire range of considered gas surface densities and fits remarkably well 
the observational data of Bigiel et al. (2008,2010) and Kennicutt (1998). This dependence 
on metallicity of the KS relation may well explain the scatter (or part of it) that is seen in 
the observationally derived relations. 
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